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Real time correlations at finite Temperature for the Ising model. 
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After having developed a method that measures real time evolution of quantum systems at a finite temperature, 
we present here the simplest field theory where this scheme can be applied to, namely the 1 + 1 Ising model. 
We will compute the probability that if a given spin is up, some other spin will be up after a time t, the whole 
system being at temperature T. We can thus study spatial correlations and relaxation times at finite T. The 
fixed points that enable the continuum real time limit can be easily found for this model. The ultimate aim is to 
get to understand real time evolution in more complicated field theories, with quantum effects such as tunneling 
at finite temperature. 
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1. Introduction 

The description of processes in real time for 
quantum systems at finite temperature is a hard 
subject that has, to my knowledge, not been fully 
understood. Most work on this field is based on 
the linear response approximation, that leads to 
computation of real time correlation amplitudes of 
operators in a "thermal state" mixture. We have 
developed a method in which one computes di- 
rectly transition probabilities among observables 
for subsets of the whole quantum system. If one 
wants to treat consistently a system at finite tem- 
perature but still measure some time correlations, 
this can only make sense for small subsystems of a 
large system in thermal equilibrium. In the ther- 
modynamic limit these correlations should con- 
verge to an answer, making the use of an external 
bath consistent with the full quantum treatment. 
We will see that this formalism gives a good prob- 
abilistic interpretation, in contrast to the one ob- 
tained from the usual treatment. 

In an earlier work j^] we have presented the for- 
malism for several problems in quantum mechan- 
ics with one or a few degrees of freedom, showing 
interesting phenomena such as real time tunnel- 
ing at finite T, but still without being able to 
show the self consistency expected in the ther- 
modynamic limit mentioned above. For this pur- 
pose and to show that the same formalism can 
be applied to field theories, the 1-1-1 dimensional 
Ising model will be computed in this work. In 



this model, also called the quantum Ising model, 
the spins in the spatial dimension will be kept 
at discrete distances while the time axis (real or 
Euclidean) will reach the continuum. In this pro- 
cess the coupling constants have to be renormal- 
ized so as to get a proper continuum limit. These 
renormalization trajectories were known for the 
Euclidean Hamiltonian formalism ||^ and will be 
obtained by analytical continuation for the real 
time case. The lack of such analytical expres- 
sion for more complicated theories, such as for 
the interesting electroweak one to study Baryon 
Number generation |3| , could turn out to be a big 
hurdle. The gauge fields can be approximated by 
discrete subgroups to get finite dimensional trans- 
fer matrices as for the spin system. 

2. Spin correlations in time 

To fix ideas let us ask for the probability, with 
the quantum Ising model, for some spin to be up 
at time t if a measured one was up at t = 0, the 
whole system being at finite T. 

In order to compute such a correlation, we just 
know that the system was in a mixed state de- 
scribed by p = exp(— Then we measure at 
t = a, spin si (without loss of generality we 
measure the first spin in a closed ring, the others 
staying undetermined), corresponding to apply- 
ing the projector (with a given si = ±1): 



P. 



|si)(si| * I 2,..JV 



(1) 
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on both sides of p. This operator evolves in t as 
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usual, U^^{Ps^pPsjUt , with Ut = exp{-iHt) . 

We measure again at a time t some spin 5*^ with 
the corresponding projectors P5. . The probabil- 
ity is then: 

Pf (si , SO = ^^Tr[Ps^ Ut {Ps, pPs, )Ut] (2) 

where we have discarded one of the Ps^ due to the 
cyclicity of the trace. The Normalization can be 
picked so that the probability for every Si to be 
up or down adds to 1. This expression can be seen 
as a proper transition probability in full analogy 
to the quantum mechanics case by going to 
the energy representation. Notice that we need 
three projectors inserted for this to happen, if one 
could commute [p, Pg J one would get the more 
usual correlations |^ which do not have a clear 
probability interpretation. In the classical limit 
both expressions would coincide. 

Equation (2) can be recast in a suitable form 
for using the path integral formalism 

± ± ± 

■ {s,..sl.\U+\..S,..Si..){..S,..Si..\Ut\si..s,..), (3) 

obtaining a product of three Green functions. 
Each of these Green functions can be calculated 
by multiplying the transfer matrix for infinites- 
imal "dt" by itself N times to obtain times t = 
(2^) * dt. This method works also for real t as we 
are considering all the contributing phases. 

The transfer matrix that takes us from one slice 
(with spins Si) to the next (spins s^), has 2" x 2" 
elements for the usual Ising model: 

M{s;s') = Y[ (4) 

i—l..n 

In order to consider this to be the short time prop- 
agator of a quantum system, one could naively 
reinterpret the couplings {Ks,Kr) as new ones 
{ArK'g , K!^ I Ar) , for which one can easily see that 
in Eq. (4) the exponent looks like an action with 
potential and kinetic energy terms, as in first 
quantization. The problem that one encounters 
in this naive approach is that this transfer ma- 
trix does not scale correctly for decreasing Ar, 



as we need to reach criticality for continuum r. 
Luckily for the asymmetric Ising model in 2 di- 
mensions one can use duality to find a critical 
line at: 'As}i\{Ks)s\v^(Kr) = 1. This was used 
to get scaling in the continuum Euclidean time 
case P, by taking Ks{At) = K'^At and then 
Kr =^arcsinh(l/ sinh(2Ar)). The new ingredi- 
ent comes in when one wants to define the trans- 
fer matrix for real times At. I assumed analytical 
behavior in the phase transition line, so that for 
Ks{At) = i K'gAt we would get criticality with 

= 1 arccosh(l/sin(2Ai)) - iir/A (5) 

obtaining a complex coupling Kt with a small 
but important phase (for scaling). With these 
couplings for the real time propagators we have 
been able to get very good convergence to the con- 
tinuum limit with a few steps of multiplying the 
transfer matrix by itself. The convergence is easy 
to test by using the composition rule for Green's 
functions, M2At = M^t ■ MAt- We still have the 
freedom to dial K'^ in order to have a stronger 
(weaker or even antiferro.) spatial binding. 

The three Green's functions needed in Eq.(3) 
can then be computed by multiplying the Eu- 
clidean or real time transfer matrices with the 
couplings as seen above. One needs N matrix 
multiplications by itself to reach times /3 = 2^Ar 
or t = 2^ At. In order to construct numerically 
the transfer matrix it is very convenient to de- 
fine first an array S(2", n) giving for each possi- 
ble spin configuration in space, the value of ev- 
ery spin. With this array Eq.(4) is easily rewrit- 
ten and Ma* (§;§') can be computed and stored 
for each pair of spin configurations s and s'. Af- 
ter multiplying the M's by themselves to obtain 
the Green's functions one still has to compute the 
product of the 3 Green's functions, making sure 
that the observables Si and Si (for a given i) have 
a fixed value and all the other ones get summed 
up. With this matrix method it is easy to get 
up to 12 or 14 spins on a spatial ring, with no 
stringent limits on t or 13 — 1/T. 

We have made various simulations in order to 
measure the probabilities, Pf{si,Si), to get a 
spin up at position i after a time t of having 
prepared si up. At high Temperatures the pre- 
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Figure 1. Probability vs. space and time to find 
the spins aligned at high T. 




Figure 2. Prob. for spin alignment at lower T. 



pared spin is initially aligned with its neighbors 
with a short spatial correlation length, the whole 
signal relaxing in time to eventually disappear 
(P = 50% everywhere) as seen in fig.l. At inter- 
mediate T's the initial spatial correlation length 
and the relaxation time are larger, as the fluc- 
tuating spins cannot screen the signal as at high 
T's. At even lower T's the spins get all almost 
aligned initially and continue in this situation for 
long times, as if going into an ordered phase, with 
interesting oscillating spin waves on top of the 
alignment, as seen in fig. 2. 

By increasing the coupling K'^, the stiffness 
grows, as if we had effectively decreased the tem- 



perature, and in fact we get similar phenomena 
as above but at already higher T's. At very weak 
K'g the prepared spin is almost uncoupled with 
the rest and mainly just oscillates down and up 
again. Changing the sign of K'^ we even get anti- 
ferromagnetic behavior among the spins. 

3. Conclusions 

We have worked out for the 1 + 1 dimensional 
quantum Ising model the spin probability corre- 
lations in real time and finite T. We find results 
that seem very plausible if one could design an 
experiment to measure these correlations. Be- 
sides the interesting results for this model, this 
work shows that this method can be applied to 
field theories as well as for simple quantum sys- 
tems. Moreover we have been able to check for 
this model that if we increase the number of spins 
on a (larger) spatial ring, we reach limiting func- 
tions for the correlations, thus reaching the ther- 
modynamic limit. In this limit the whole scheme 
is consistent in the sense that the initial prepara- 
tion of spin si does not take the full system out 
of thermodynamic equilibrium and the rest of the 
system can be considered as a genuine quantum 
heat bath to study the time evolution. 

Our ultimate goal is to find a way to treat the 
field theory case of tunneling with instantons in 
real time and finite T. As we have seen, it is in 
principle doable with these methods but will be 
hard to implement due to the number of degrees 
of freedom after discretizing the gauge fields. The 
other hard part will consist in finding the renor- 
malization trajectories for the coupling, which 
could lie in the compex plane, in order to reach a 
sensical continuum time limit. 

I would like to thank J. Polonyi, P. Rujan and 
L. PoUey for useful discussions on this work. 
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